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We present a first-principles study of electronic and magnetic properties of thin Co films on a 
BaTiO3(001) single crystal. The crystalline structure of 1, 2, and 3 monolayer thick Co hlms was 
determined and served as input for calculations of the electronic and magnetic properties of the 
hlms. The estimation of exchange constants indicates that the Co hlms are ferromagnetic with 
a high critical temperature, which depends on the him thickness and the interface geometry. In 
addition, we calculated x-ray absorption spectra, related magnetic circular dichroism (XMCD) and 
linear dichroism (XLD) of the Co L 2,3 edges as a function of Co him thickness and ferroelectric 
polarization of BaTi 03 . We found characteristic features, which depend strongly on the magnetic 
properties and the structure of the him. While there is only a weak dependence of XMCD spectra 
on the ferroelectric polarization, the XLD of the hlms is much more sensitive to the polarization 
switching, which could possibly be observed experimentally. 

PACS numbers: 75.50.Cc, 71.20.Lp, 71.15.Rf 


The interface between a magnetic thin him and a fer¬ 
roelectric material is the subject of several recent in¬ 
vestigations. The major part of these studies is de¬ 
voted to Fe/ferroelectric interfaces,since 
ferromagnetic iron is supposed to be a good candidate 
as a ferromagnetic electrode in two-component multi- 
ferroics. Although it was shown, that a ferroelectric 
him can be grown on an iron substrate and a func¬ 
tional heterojunction can be fabricated,the multi- 
ferroic ehFects by polarization switching were found to be 
not strong enough, since the main changes of the func¬ 
tional properties occur only in the vicinity of the inter¬ 
face.Another impediment is non-ideal compatibility 
between an iron and an oxide surface. Until now it was 
not shown that junctions with symmetric interfaces, a 
highly desirable condition for coherent electronic trans¬ 
port,^*^ can be grown. Accordingly, ferromagnetic oxides 
Lai-a^Sr^^MnOs and SrRuOs were used as electrodes in 
multiferroic junctions. Such interfaces can be al¬ 
most ideally grown, but the Gurie temperature of these 
oxides is too low for functional devices. Therefore, search 
for an appropriate ferromagnetic/ferroelectric interface is 
still in progress. 

Among the above cited studies only few investigations 
deal with ultrathin ferromagnetic hlms on a ferroelec¬ 
tric substrate. 16 particular, it was shown that 

two monolayer thick Fe hlms on the BaTiO3(001) surface 
might not be ferromagnetic because of the him geometry 
and magneto-elastic properties of iron.^ Additionally, a 
substantial charge and spin moment transfer was found 
at the interface by altering the polarization direction.^ 

Thus, ultrathin hlms of Fe on a BaTiO3(001) sin¬ 


gle crystal are magnetically unstable,^ but cobalt ex¬ 
hibits usually stable ferromagnetic characteristics in 
many nanostructures. Therefore, we continue our work 
on ultrathin metallic hlms on ferroelectric single crystals 
and suggest in this paper to use Go as the ferromagnetic 
material on BaTiO 3 ( 001 ). 

Despite the comprehensive review about the progress 
in this held given by Vaz et al?^ there are only few 
information about the interface of Co and perovskites 
like BaTi 03 (BTO). In the framework of spin-polarized 
DFT calculations, as implemented in the Vienna ab ini¬ 
tio Simulation Package (VASP), multiferroic tunnel junc¬ 
tions of Co/BTO/Co were investigated by Gao et al}^ 
They showed that a critical thickness of BTO unit cells 
is necessary for the appearance of ferroelectricity, which 
is inhibited by a depolarizing electrostatic held, caused 
by dipoles at the ferroelectric-metal interfaces.In the 
work of Oleinik et al?^ hrst-principles calculations were 
applied on Co/STO/Co(001) magnetic tunnel junctions, 
where a strong covalent bond between Co and O and 
an induced magnetic moment of |ms| = 0.25/ip at the 
Ti atom was observed. This is similar to the case of 
Fe/BTO/Fe tunnel junctions. In another work, Polisetty 
et al?^ applied piezoelectrically controlled strain for elec¬ 
tric tuning of exchange-bias helds of BTO/Co/CoO het¬ 
erostructures. 

In our hrst-principles study, we investigated systemat¬ 
ically the crystalline structure of ultrathin Co hlms on 
BTO and their electronic and magnetic properties in de¬ 
pendence on the him thickness and the polarization of 
BTO. The calculations were performed using a so-called 
multi-code approach, in which atomic positions were ob- 
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tained using a pseudo-potential code, VASP.^^’^^ This 
information serves as input for electronic and magnetic 
structure calculations with different multiple-scattering 
Green function methods. 

In addition, we simulated x-ray absorption spec¬ 
tra (XAS) and the related x-ray magnetic circular dichro- 
ism (XMCD), which is the method of choice to prove ex¬ 
perimentally the change of the magnetic structure. The 
x-ray linear dichroism (XLD) is dicussed in respect to 
the occupation of d orbitals. These methods are local 
and site-sensitive and offer the opportunity to investi¬ 
gate both the magnetic and structural properties as a 
function of layer thickness.So, we traced the depen¬ 
dence of XAS of Co L 2,3 edges as a function of the Co 
film thickness and the polarization direction of BTO. All 
results were compared with the previous results for Fe 
thin films on BTO.^’^^ 

Our paper is organized as follows. Our multi-code ap¬ 
proach and computational details are presented in the 
next section I. Then, in sections II and III, we discuss 
structural, electronic, and magnetic properties of ultra- 
thin Co films on BTO. Sections IV and V deal with com¬ 
putational simulations of XAS, XMCD, and XLD spec¬ 
tra. Conclusion is offered in section VI. 


I. COMPUTATIONAL DETAILS 
A. Structural optimization 

The information about the crystalline structure of 
Co/BTO(001) was obtained using projector-augmented 
wave pseudo-potentials^^ implemented in the VASP 
code.^^’^^ The plane-wave basis was taken with a cut¬ 
off energy of 400 eV. The calculations were performed 
within the local spin-density approximation (LSDA). We 
used the parametrization of Perdew and Wang for the 
exchange-correlation functional (PW-functional).^^ Here, 
we believe that the impact of electronic correlations on 
the interface magnetoelectric coupling is minor. 

We performed the structural optimization of 
Co/BTO(001) in dependence on the Co film thick¬ 
ness and the ferroelectric polarization P of BTO (001). 
To model the (001) surface of polar BTO, we constructed 
a 5 unit cells (^ 2nm) thick supercell, with a vacuum 
spacer of 2nm along the [001] direction. The lattice 
parameters were set to the theoretical equilibrium values 
of tetragonal BTO:^ a = 3.943 A and c/a = 1.013. The 
intralayer cation-oxygen displacements 5 = zq — ^^cation 
(see Fig. 1) in BTO cause the ferroelectric polarization, 
P, along [001]. When P is antiparallel to the surface 
normal (P|), the oxygen in each monolayer (ML) are 
higher then the cations {6 > 0), and vice versa, the 
ferroelectric state means < 0. Since the (001) 
surface of BTO is Ti02 terminated^ for both directions 
of P, the Co atoms find their relaxed positions atop 
the oxygen atoms. This is in agreement with the recent 
experimental data for Fe/BTO(001).^^ Each Co mono¬ 


layer of the slab contains, therefore, two atoms per unit 
cell while the Co film thickness, L, varies between one 
and three ML. For Co/BTO(001) the atomic positions 
of the four top BTO layers and all Co ML were allowed 
to relax. After the relaxation, the calculated forces 
are always less than 0.5 x 10“^eV/A. The Brillouin 
zone of the slab was sampled with a 10x10x6 /c-point 
Monkhorst-Pack mesh^^ during the force minimization. 

B. Electronic and magnetic structure calculations 

The relaxed geometry was used for further first- 
principles calculations using the multiple-scattering 
Green function method (Hutsepot) within the atomic 
sphere approximation to the crystal potential.We 
took an angular momentum cutoff of /max = 3 for the 
Green function, a A;-point mesh of 24 x 24 x 12 for the BZ 
integration and 24 Gaussian quadrature points for com¬ 
plex energy contour integration. With the Green function 
G{E) of the system, all quantities of interest follow in a 
straightforward way. In particular, we are interested in 
ground state properties like the density of states (DOS) 
and the local magnetic moments^^ in dependence on Co 
film thickness L and ferroelectric polarization direction 
P of BTO. 

To describe magnetic properties of Co/BTO(001), we 
calculated interatomic exchange coupling parameters Jij 
using the magnetic force theorem implemented within the 
Green function method.The exchange coupling con¬ 
stants Jij can be used to obtain spin-wave spectra by the 
diagonalization of the Heisenberg Hamiltonian 

H = -'"^JijSiej - ^ Pmae (ef)^ , (1) 

i/j i 

where i and j label magnetic atoms, is a unit vec¬ 
tor in the direction of the magnetic moment of the i- 
th atom and Pmae is the magnetic anisotropy energy 
(MAE), which is positive (negative) for the case of the 
easy-axis (easy-plane) anisotropy type. It can be calcu¬ 
lated directly with the KKR method by using again the 
magnetic force theorem.The critical temperatures were 
estimated using a Monte Carlo (MC) simulation with the 
model Hamiltonian in Eq. (1). Eor the simulation, two 
dimensional supercells, which repeat the unit cell (see 
Eig. 1) 60 X 60, 80 x 80, and 100 x 100 times in x and y 
direction were constructed. In those directions, we con¬ 
sidered also periodic boundary conditions and restricted 
the calculation only to the magnetic atoms of the unit 
cell (2 Co atoms for L = 1, 4 Co atoms for L = 2, and 6 
Co atoms for L = 3). Within the supercell, the magnetic 
moment at lattice site i interacts with its neighbors at 
the site j via the first-principles Jij. During a MC run, 
a lattice site j with the magnetic moment vector Cj was 
chosen and a new random direction e' was created. The 
energy of the system determines whether e' or Cj was 
kept. Performing this procedure N times on a lattice of 
N sites is defined as one MC step. The starting point 
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FIG. 1: Geometry of Go films on the BaTiOa surface for 1 ML 
(a), 2 ML (b), and 3 ML (c) cobalt thicknesses, respectively. 
Spheres represent Go (brown), Ba (gray), Ti (blue), and oxy¬ 
gen (red) sites. Values of selected interatomic distances for 
polarizations P = Pi,Pt presented in Table 1. The pic¬ 
ture of the underlying structure was done with VESTA.The 
numbers for the Go atoms symbolize the atoms of the unit cell 
and correspond to those in the figures of the supplementary 
material. 


was a high-temperature disordered state above the criti¬ 
cal temperature Tq. Iu the course of the simulatious, the 
temperature was stepwise reduced uutil maguetic order- 
iug was reached. For each temperature T, the thermal 
equilibrium was assumed to be reached after 20 000 MC 
steps. The thermal averages were determiued over 20 000 
additioual MC steps. Tq was theu obtaiued from the 
temperature depeudeucy of the maguetic susceptibility. 
With respect to fittiug procedures aud fiuite temperature 
sampliug, all three supercells yielded the same critical 
temperature withiu au uucertaiuty rauge of ±5K. More 
details of our MC scheme cau be fouud iu Refs. 38-40. 


C. X-ray absorption spectra simulations 

Usiug the self-cousisteut potentials of the scalar- 
relativistic Greeu fuuctiou method (Hutsepot), the cal- 
culatious were extended to a spin-resolved fully rela¬ 
tivistic KKR code (SPR-KKR).^^ This Green function 
method^^ allows the calculation of x-ray absorption co¬ 
efficients iJi^{E) in dependence on the energy E and the 
polarization A of the x-rays.For SPR-KKR calcu¬ 
lations of magnetic moments, XAS and DOS, we have 
taken 144, 578 and 225 /c-points in the irreducible part 
of the Brillouin zone, respectively. The calculation of the 
structure constants were done using the Ewald method 
with parameter 2.0 and 4.0 in real space and reciprocal 
space, respectively. For the Ewald parameter connecting 
the summation in real and reciprocal space a value of 
0.8 was used, making the summation in the real space 


TABLE I: Relaxed interatomic distances (in A) of 
(Go)l/BTO( 001) (L < 3 and P = For each L, the 

distances between Go ML i and j are denoted by dij while 
the ^-separation from interfacial O (Ti) to Go of the first 
(second) ML is shown by doi {dT2)- For comparison, also 
nearest neighbor Go distances in bulk (hep) from Ref. 45 and 
distances for (Fe)L/BTO(001) are shown. 



Co 

Pi Pi 

Fe 

Pi Pi 

L = 1 doi 

1.778 

1.784 

1.774 

1.781 

II 

to 

0 

1.856 

1.853 

1.857 

1.855 

di2 

1.116 

1.107 

1.049 

1.054 

dT2 

3.022 

2.927 

2.971 

2.918 

0 

CO 

II 

1.827 

1.829 

1.843 

1.849 

di2 

1.183 

1.179 

1.218 

1.241 

d23 

1.154 

1.153 

1.134 

1.114 

dT2 

3.094 

3.014 

3.214 

3.121 

bulk (hep): 

2.478 




converge faster. These parameters were applied for all 
considered Co/BTO systems. 


II. CRYSTALLINE STRUCTURE OF 
ULTRATHIN CO FILMS ON BTO 

For the ultrathin Co films deposited on BTO(OOl), 
the optimized interatomic distances are collected in 
Tab. I, in comparison with the corresponding structure 
of Fe/BTO(001).^’^^ The definition of our distances are 
given in Fig. 1. Below in the text, the interface Ti02 
layer is denoted as I for all films. The Co layers towards 
the surface are labeled I + 1, 1 + 2, and 1 + 3. 

Structural relaxation, with the in-plane degrees of free¬ 
dom for the layers starting with / —2, energetically favors 
the structure with the Co sites in layer I + 1 on top of 
the O sites in layer 1. In the case of 1 ML Co (L = 1), 
the structure is similar to that of Fe/BTO(001).^’^^ For 
both multiferroic systems, the polarization state P^ re¬ 
sults in a slightly shorter distance, doi, between Co (Fe) 
and interfacial O. 

When the Co film thickness increases to L = 2 and 3, a 
distorted bodycentered tetragonal (bet) lattice is formed, 
although it is not typical for Co.^^ We found that the po¬ 
larization reversal affects mainly the positions of interfa¬ 
cial Ti (I) and consequently those of the Co (1 + 2) atoms, 
labeled as Copa and Coxi, respectively (see Fig. 1). This 
is important concerning the magnetic properties of the 
films. 


III. ELECTRONIC AND MAGNETIC 
PROPERTIES OF ULTRATHIN CO FILMS ON 
BTO 

In order to determine the magnetic ground state for 
the three film thicknesses, we calculated the total en- 
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TABLE II: SPR-KKR results of magnetic spin moments ttIq 
(/ie/atom) for (Co)l/BTO( 001) in dependence on the Co film 
thickness L and the direction of ferroelectric polarization P 
{M II [001]). The Co labels follow from Fig. 1. In I + 1 and 
1 + 3 the Co atoms are equivalent. 




L = 

n 

1 

Pt 

L = 
Pi 

2 

Pt 

L = 

Pi 

3 

Pt 

Co 

1 + 3 





1.39 

1.42 

CoBa 

1 + 2 



2.00 

2.03 

1.60 

1.64 

Coxi 

1 + 2 



1.89 

1.93 

1.45 

1.47 

Co 

I + l 

1.71 

1.71 

1.62 

1.62 

1.42 

1.37 

0 

I 

0.08 

0.08 

0.03 

0.03 

0.06 

0.05 

Ti 

I 

-0.11 

-0.18 

-0.12 

-0.16 

-0.03 

-0.08 


Co bulk (hep): 1.62 


TABLE III: SPR-KKR results of magnetic orbital moments 
Trio (//B/atom) for (Co)l/BTO( 001) in dependence on the Co 
film thickness L and the direction of ferroelectric polarization 
P (AT II [001]). The Co labels follow from Fig. 1. In I + 1 
and 1 + 3 the Co atoms are equivalent. 




L = 

Pi 

1 

Pt 

L = 

Pi 

2 

Pt 

L = 

Pi 

3 

Pt 

Co 

1 + 3 





0.09 

0.08 

COBa 

1 + 2 



0.08 

0.08 

0.07 

0.07 

Coxi 

1 + 2 



0.09 

0.08 

0.07 

0.07 

Co 

I+l 

0.17 

0.17 

0.07 

0.06 

0.09 

0.09 

Co bulk (hep): 

0.06 







ergy for the ferromagnetic (EM) and antiferromagnetic 
(AFM) configurations of the magnetic moments of the Co 
atoms. In all cases, the EM solution was preferable, inde¬ 
pendent of thickness or polarization. So, the thin cobalt 
films are strong ferromagnets and show a lower sensi¬ 
tivity to structural transformations then the iron layers, 
which showed a considerable change in the magnetic or¬ 
der for different layer thicknesses. For the 2 ML Fe film 
on BTO, the mpe in layer I+l was almost quenched while 
the sizable moments in the surface layer 1 + 2 are ordered 
antiparallel. This results in a total magnetic moment of 
m ^ 0. Deposition of a third Fe monolayer restored the 
ferromagnetic order of the 1 ML Fe film. 

A. Magnetic moments and density of states 

The calculated spin moments (see Tab. II) for the Co 
atoms show a strong dependence on the geometry of the 
films and the hybridization of Co 3d states with the elec¬ 
tronic bands of the substrate. The different environment 
of e.g. CoBa and Coxi, in particular, the atomic vol¬ 
umes and the band hybridization, influences also these 
moments. 

In the case of L = 1 and L = 2, the spin moments 
for Co atoms at the surface, are larger then for the Co 
bulk (hep) because of the symmetry reduction and an 
enhanced volume per Co atom. The Co magnetic mo- 


TABLE IV: SPR-KKR results of magnetic orbital moments 
rrio (/xe/atom) with M || [100]. The Co labels follow from 
Fig. 1. In this held direction the Co atoms are nonequivalent. 




L = 

Pi 

1 

Pt 

L = 

Pi 

2 

Pt 

L = 

Pi 

3 

Pt 

C06 

1 + 3 





0.12 

0.07 

C05 

1 + 3 





0.07 

0.11 

C04 

1 + 2 



0.12 

0.12 

0.07 

0.07 

C03 

1 + 2 



0.12 

0.13 

0.06 

0.06 

C02 

I+l 

0.16 

0.12 

0.08 

0.08 

0.06 

0.04 

Coi 

I+l 

0.13 

0.17 

0.07 

0.06 

0.07 

0.05 



E-E^ [eV] E-E^ [eV] E-E^ [eV] 


FIG. 2: Spin-resolved (t, i) density of d states for Co and Ti 
atoms in dependence on him thickness L. It is only shown for 
Pf, because there are only small diherences for P^. 


ments in layer I + 1 for L = 2 approach their bulk value. 
For L = 3, the spin moments of Co (I + 1) and Coxi are 
smaller in comparison to L = 2, because of slightly elon¬ 
gated distances di 2 and dx 2 : respectively. In the case of 
L = 3, we observe for all three layers a quenching of the 
spin moment in comparison to Co bulk (hep), which re¬ 
sults from a strong reduction of the volume per Co atom. 

Besides, we also observe induced magnetic moments at 
the Ti and O atoms at the interface (see Tab. II). For the 
Ti atoms, these induced magnetic moments are antipar¬ 
allel oriented to the direction of the magnetic moments 
of the Co atoms and are originated, similar to the case 
of Fe/BTO(001), from the hybridization of Ti and Co 
3d states (see Fig. 2). The strong interaction between 
the host and the Co films widens the DOS of the layers 
at the interface, while the 3d states of Co become very 
narrow towards the surface, which results from the re¬ 
duced coordination number of Co atoms at the surface 
(in comparison to Co bulk). The impact of the crystalline 
environment is very strong on the DOS and magnetic 
moments of the Co atoms. As the result, the magnetic 
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Layer: I + 1 



^Ba 

^T, 


a o 


Layer: I + 1 


Layer: 1 + 2 



^ 100 ] 


FIG. 3: Calculated intralayer exchange constants (in meV) 
for the I + 1 layer in 1 ML Co/BTO(001) and Pf, viewed 
from the z direction. Atoms in lighter colors are below those 
in darker colors. The numbers of the Co atoms correspond to 
those used in Fig. 1(a). 


moments on Coxi are smaller by 0.08 /tb to 0.15 /tb then 
on CoBa(I + 2) due to the strong interaction between 
the Coxi and the Ti states. A short distance between 
Co(I + 3) and Co(I + 2) and a strong Co-Co hybridiza¬ 
tion in the case of L = 3 leads to a substantial reduction 
of magnetic moments of the surface Co layer. 

We have calculated the orbital moment of Co atoms 
along M II [001] (see Tab. Ill) and along M || [100] (see 
Tab. IV). In the last case, the two Co atoms in each layer 
are nonequivalent. In the surface layer the atoms show 
a strong dependence on the electric polarization of the 
substrate. For M || [001] the influence of the polarization 
is weak and is as large as for the Fe/BTO system. 

Duan et predicted that the magnetic anisotropy 
of thin magnetic films may be effected by the polariza¬ 
tion of the ferroelectric substrate. His calculations done 
for 1 ML Fe/BTO gave differences for the orbital mag¬ 
netic moment Arrio = mo[001] — mo[100] of 0.035 /ip and 
0.021 /ip for P| and P^, respectively. In our former work 
of 1 ML Fe/BTO, we found also an averaged anisotropy of 
the orbital magnetic moment Arrio of 0.04 /ip and 0.06 /ip 
(note that there are two Fe atoms in the unit cell of 
1 ML Fe).^^ In case of 1 ML Co/BTO, we get the av¬ 
eraged Arrio = 0.04/ip for the polarization which is 
similar to 1 ML Fe/BTO. There was no strong depen¬ 
dence on the polarization on the BTO substrate in our 
calculations. The anisotropies for the orbital moment will 
influence the preference of the easy axis of magnetization 
as discussed in the next section. 

On the other hand, for thicker Co films, the electronic 
structure can be strongly affected by a Schottky barrier, 
which develops at a semiconductor-metal interface and 
can lead to a spurious charge transfer across the inter¬ 
face.^^ However, in very thin metallic films on a semicon¬ 
ductor, the Schottky barrier is not significant and can not 
substantially modify their electronic and magnetic struc¬ 
ture. 


Q Ba o Ti OO QCo OCoTi 

FIG. 4: Calculated exchange constants (in meV) for 2 ML 
Co/BTO(001) and P^. The notation is the same as in Fig. 3 
and the numbers of the Co atoms correspond to those used in 
Fig. 1(b). Intralayer interactions in the I + 1 and 1 + 2 layers 
are separated from the interlayer interactions shown in the 
side view. The dashed lines in the two upper panels indicate 
the the plane which is shown in the side view. For symmetric 
interactions only one number is given. 


B. Magnetic interactions and Curie temperatures 

To study magnetic interactions in ultrathin Co films on 
BTO(OOl), we computed the exchange constants which 
enter the Heisenberg model (1). The most significant 
magnetic exchange interactions for Co/BTO (001) are 
shown in Figs. 3 to 5 for P = 1, 2, 3 ML, respectively, 
while for the sake of completeness all calculated coupling 
constants are presented in the supplementary material. 
They are always separated in interlayer and intralayer 
contributions. 

The key feature of the elaborated results is a very 
strong magnetic coupling between the nearest neighbors. 
This means for the 1 ML Co film the nearest neigh¬ 
bor coupling is in-plane (nearest neighbor Co distance 
is 2.79 A) with about 21meV (see Fig. 3). In this case 
the electronic density is mainly distributed within the 
x^-plane and the exchange integrals represent an overlap 
between the electronic wave functions. Hence, compar¬ 
ing with Co bulk (hep) the magnetic coupling in 1 ML 
remains stronger despite the larger nearest neighbor dis¬ 
tance (in Co bulk (hep) the nearest neighbor magnetic 
coupling is ~ 13meV). 
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Layer: I + 1 


Layer: I + 2 


Layer: I + 3 




Side view 



^Ba 


O O 

Q Co 
^ COTi 


FIG. 5: Calculated exchange constants (in meV) for 3 ML Co/BTO(001) and P^. The notation is the same as in Fig. 3 and 
the numbers of the Co atoms correspond to those used in Fig. 1(b). Intralayer interactions in the three layers I + 1, I + 2 and 
1 + 3 are separated from the interlayer interactions shown in the side view. The dashed lines in three upper panels indicate the 
the plane which is shown in the side view. For symmetric interactions only one number is written (see text). 


In contrast, for the films with an out-of-plane com¬ 
ponent (2 ML, 3 ML), the distances towards the Co 
atoms in the next layer are shorter (2.23 A). Therefore, 
the nearest neighbor coupling is now between those lay¬ 
ers and the wave function overlap is even stronger. It 
leads to a strong magnetic interaction between 32meV 
to 42meV for 2 ML (see Fig. 4). In the case of 3 ML, 
the electronic density redistributes in all directions and, 
therefore, the strength of the interlayer exchange interac¬ 
tions with respect to those of 2 ML is reduced to 13 meV 
to 20meV (see Fig. 5). Comparing these results with 
the exchange constants in Co bulk, we found in general 
a strong asymmetry between in-plane and out-of-plane 
components. The latter are larger while the in-plane di¬ 
rections are strongly reduced (see Figs. 4 and 5). 

Furthermore, we obtain a partial mediation of the mag¬ 
netic coupling by the BTO host. The exchange interac¬ 
tion between e.g. the second nearest neighbor Co atoms 
depends on the underlying atom, either ^ —3.3 meV or 
^ 0.4 meV with mediating Ba or Ti, respectively (see 
Fig. 3). This fact can also be observed for 2 and 3 ML 
(see Figs. 4 and 5) and is evident from the comparison 
of the calculated Jij for supported and unsupported thin 
Co films, shown in the supplementary material. While 
using the same geometry for the Co films, we removed 
the BTO substrate in the calculations. Immediately, the 


TABLE V: Magnetic anisotropy energy and critical tempera¬ 
ture for (Co) L /BTO(OOl) in dependence on the Co film thick¬ 
ness L and the direction of ferroelectric polarization P. Tc 
in brackets is for Emae = 0 meV. 


L 

-Emae 

A 

(meV) 

A 

A(K) 

A A 

1 

-1.847 

-1.483 

360 (307) 

298 (250) 

2 

0.680 

0.958 

844 (820) 

818 (780) 

3 

-0.270 

-0.438 

567 (550) 

580 (550) 


symmetry of x and y directions returns and some ex¬ 
change interaction values differ more than 10 meV. 

For the Monte Carlo simulations with the classical 
Heisenberg Hamiltonian (1), we have also computed the 
magneto-crystalline anisotropy energy (MAE) E^mae = 
E^[ooi] — E^[ioo] foi* various Co thicknesses and polariza¬ 
tions (see Table V). We found an out-of-plane magnetiza¬ 
tion direction for the cases of 1 and 3 ML Co/BTO(001), 
while in the case of 2 ML Co/BTO(001), the magnetiza¬ 
tion direction is in-plane. The change in the MAE with 
the polarization switch can reach up to 25 % to 30 %. It 
is remarkable, that in the 2 ML Co case, the magnetiza¬ 
tion direction is lying within the x^-plane, which can be 
explained by the small interlayer distance (see Table I). 

Together with the calculated Jij parameters, those 
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FIG. 6: Calculated XAS (upper and middle part) and 

XMCD (lower part) of 1 ML Co on BTO(OOl) (P^) in com¬ 
parison to Co bulk (hep). Open arrows indicate difference in 
peak position (see text). 


MAE values were used to determine the critical tem¬ 
peratures Tq in dependence on the layer thickness L 
and polarization P (see Table V). In general, the thin 
(Co)i,/BTO films are ferromagnetic at room tempera¬ 
ture for 1 / > 2 . 

The Tq increases not monotonically with the thickness, 
as expected from the large change of the magnetic cou¬ 
pling parameters. The main contribution originates from 
the strong interlayer coupling. The different polarization 
directions have only a small influence to the the value 
of Tq which follows from the similar Jij parameters for 
or P| (for all values see the supplementary material). 
As expected, a value for MAE |Pmae| > 0 increases the 
Curie temperature (see Table V). 


IV. XAS AND XMCD 

The Co L 2,3 edges are in the focus of our XMCD dis¬ 
cussion. We have calculated x-ray absorption spectra of 
right- and left-circular polarized x-ray radiation in de¬ 
pendence on the Co film thickness L and the electric 
polarization direction P of BTO. The difference of the 
absorption coefficients A/i = /i+(P)—/i“(P) = XMCD is 
normalized to the number of Co atoms in the multilayers. 
All x-ray absorption spectra were broadened by Lorentz 


FIG. 7: Calculated XMCD of Co L 2,3 edges in Co/BTO in 
dependence on film thickness and electric polarization P of 
BTO {P^ solid lines, P^ dotted lines). Arrows indicate differ¬ 
ent small features (see text). The thin vertical line visualizes 
the difference in the peak positions. 


convolution with a core hole width of 0.9 eV and 0.25 eV 
for the L 2 and L 3 edge, respectively. In our simulations 
of XAS and XMCD the magnetization AT, the electric 
polarization P and the incidence of light are parallel to 
the z-axis (surface normal) of the Co/BTO(001) system. 

In Fig. 6 the calculated x-ray spectra /i^(E) (upper 
part) and the related difference spectra (lower part) are 
shown for 1 ML Co/BTO(001) (P^) in comparison to 
calculations performed for Co bulk (hep). The results of 
our XMCD calculations for (Co)l/BTO with L = 1,2,3 
and P^ and P^ are summarized in Fig. 7. 

In all cases, we observe the well-known energy depen¬ 
dence of XMCD (Fig. 7), where the intensity at the 
L 3 (L 2 ) edge is decreased (enhanced) due to different se¬ 
lection rules of spin-up and spin-down electrons in the 
ferromagnetic phase. The behavior is opposite at the L 3 
(L 2 ) edge in case of /i“. The dependence on the elec¬ 
tric polarization P of BTO is in all cases very weak (see 
dotted lines in Fig. 7). 

In comparison to Co bulk, we found small but signif¬ 
icant differences concerning ultrathin films. Note that 
the maxima of the L edges of left and right circularly 
polarized light are not at the same energetic position as 
observed for Co bulk XAS (see open arrows in Fig. 6 ). 
This can be attributed (as in the case of Fe/BTO^^) to 
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FIG. 8: Calculated XLD of Co L 2,3 edges in Co/BTO for 
polarization P (P^ solid lines, Pi dashed lines). The thin 
vertical line visualizes the difference in the peak positions. 


the cancellation of degenerate d states at the interface 
or/and surface of the film. 

In the energy dependence of the L 3 edge, we observe a 
small structure of the XMCD peak at higher energy (see 
closed arrows in Fig. 7) in dependence on layer thick¬ 
ness. The appearance of this structure can also be ex¬ 
plained if we consider the spin- and layer-resolved DOS 
(not shown here) of d states as demonstrated for the 
system Fe/BTO.^^ In general, it is possible to separate 
all the different contributions by means of layer-resolved 
XMCD calculations. This dependence on layer thickness 
should show up in corresponding experiments. 


V. X(M)LD 

X-ray absorption linear dichroism can be applied in 
different modes.In nonmagnetic materials with cubic 
symmetry the x-ray absorption intensity is independent 
on the orientation of the polarization vector of light {u) 
relative to the sample. As the symmetry is lowered (like 
in a film), the intensity is directly proportional to the 
number of empty valence states in the direction of u. The 
polarization vector u acts like a search light for the direc¬ 
tion of maximum and minimum empty valence states. 

In case of a magnetic sample with cubic symmetry, the 
spin-orbit coupling leads to magnetic linear dichroism. 


where the x-ray absorption intensity is different for u 
aligned parallel and perpendicular to the spin direction. 
This kind of dichroism is often considered in the discus¬ 
sion of magneto crystalline effects (for example in Ref. 
49). 

In our magnetic ultrathin Co films on BTO we have 
both natural and magnetic linear dichroism effects. First, 
we present results for a fixed direction of magnetization 
M II [001] and linearly polarized light with polarization 
along [001] (/i^) and [100] (/i^^), respectively. The differ¬ 
ence XLD = — iP{E) between both x-ray absorp¬ 

tion spectra is shown in Fig. 8 . 

The strong dependence of the XLD on film thickness 
can be related to the occupation of d orbitals. A de¬ 
tailed analysis is possible concerning the single (pj^rrij) 
contributions of the initial P 3/2 {Pi/ 2 ) state to the p ^ d 
transitions of the different Co absorbers in the film. In 
the 1 ML case, we get the main positive contribution 
from the transition of rrij = — 1/2 into d^z‘^_r 2 states. 
The contributions of dxz and dyz are small because of 
the strong interaction between Co and Ti atoms at the 
interface. The XLD of 2 ML reflects the change of occu¬ 
pancy between in-plane and out-of-plane orbitals. In the 
3 ML case, we found that the behavior is similar to the 
Co bulk, but with much more fine structures. These fine 
structures are related to the lifting of degeneracy in the 
final d state in the thin film. 

Besides, we found that the XLD is much more sensi¬ 
tive concerning the dependence on the polarization of the 
BTO substrate as the XMCD. This emphasize how the 
reversal of polarization P changes the occupation of d 
states. 

We have also investigated x-ray absorption in depen¬ 
dence on the magnetization direction M for fixed polar¬ 
ization u II [ 001 ] but spin-orbit interaction and crystal 
field are both to weak to eliminate magneto crystalline 
anisotropy. 


VI. CONCLUSION 

In this first-principles study of structural, electronic, 
and magnetic properties of ultrathin Co layers on the 
BaTiO 3 ( 001 ) substrate, we showed, that the crystalline 
structure of (Co)l/BTO(001) {L = 1,2,3) interfaces is 
very similar to the one of (Fe)L/BTO(001) films pre¬ 
sented in our previous study. ^ Additionally, we found also 
for the Co layers only a small dependence on the polariza¬ 
tion direction of the substrate for all investigated prop¬ 
erties. However, in contrary to the Fe/BTO(001) case, 
Co films on BTO (001) exhibit a stable ferromagnetic or¬ 
dering at room temperature which depends strongly on 
the layer thickness and is the strongest for L = 2. The 
magnetic interaction is mainly featured by the coupling 
between the nearest neighbors. For L > 2 cases, the 
strongest interactions arise between the adjacent layers, 
while the intralayer magnetic coupling was found to be 
rather weak. This results from short interlayer distances 
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leading to a strong hybridization between Co 3d electrons 
of adjacent layers. Surprisingly, the easy axis turns for 
2 ML from a out-of-plane magnetization (in L = 1 and 3) 
to a in-plane magnetization, which was obtained from the 
direct calculation of the MAE or the calculated orbital 
moments. 

For a comparison to possible experimental measure¬ 
ments, we simulated x-ray absorption spectra and related 
x-ray magnetic circular and (magnetic) linear dichroism 
to trace the change of the spectra under the polariza¬ 
tion switching in the BTO substrate. While the XMCD 
depends only weakly on the substrate polarization, sim¬ 
ilar to our previous study for (Fe)L/BTO(001) films, 
the XLD shows indeed significant changes under polar¬ 


ization switching, which can probably be observed exper¬ 
imentally. 
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